
PCT/GB 2004 / 0 0 0 1 2 6 




INVESTOR IN PEOPLE 




The Patent Office 



Concept House 



PRIORITY 
DOCUMENT 



Cardiff Road 
Newport 



SUBMITTED OR TRANSMITTED IN 
COMPLIANCE WITH RULE 17.1(a) OR (b) 



South Wales 
NP10 8 QO 



recd 2 3 FEB 2004 



WIPO 



PCT 



I, the undersigned, being an officer duly authorised in accordance with Section 74(1) and (4) 
of the Deregulation & Contracting Out Act 1994, to sign and issue certificates on behalf of the 
Comptroller-General, hereby certify that annexed hereto is a true copy of the documents as 
originally filed in connection with the patent application identified therein. 



In accordance with the Patents (Companies Re-registration) Rules 1982, if a company named 
in this certificate and any accompanying documents has re-registered under the Companies Act 
1980 with the same name as that with which it was registered immediately before re- 
registration save for the substitution as, or inclusion as, the last part of the name of the words 
"public limited company" or their equivalents in Welsh, references to the name of the company 
in this certificate and any accompanying documents shall be treated as references to the name 
with which it is so re-registered. 



In accordance with the rules, the words "public limited company" may be replaced by p. I.e. , 
pic, P.L.C. or PLC. 



Re-registration under the Companies Act does not constitute a new legal entity but merely 
subjects the company to certain additional company law rules. 




An Executive Agency of the Department of Trade and Industry 



BEST AVAILABLE COPY 



Dated 13 February 2004 



Signed 




* 

pole 16) 



Form 1/77 

977 



THE PATENT OFFICE 
F 



FEB 201 
NEWPORT 




.ho 



Request for grant of a patent 

(See the notes on the back of tins form. Yon can also get an 
explanatory leaflet from the Patent Office to help you SB in 
this form) 



The Patent Office 

CardifFRoad 
Newport 
South Wales 
NP9 1RH 



1. Your reference 


2QgJ WIJN/67977/000 




2. Patent application number 

(The Patent Office will Gil in this part) 


0302898.2 


P01/7700 0.00-030289B.2 



each applicant (underline all surnames) 



ABB Offshore Systems Limited 

2 High Street 

Nailsea 

Bristol BS481BS 



Patents ADP number (if you know it) ' CF=f fe^ 9^*3- I OO ( 

If the applicant is a corporate body, give the United Kingdom (GB) 

country/state of its incorporation 



4. Title of the invention 



ESTIMATING THE TIME OF ARRIVAL OF A SEISMIC 
WAVE 



5. Name of your agent (if you have one) PAGE HARG RAVE 

« AJS * . „ . Tr Southgate, Whitefriars 

Address for service in the United Kingdom Lewins Mead 

to which all correspondence should be sent BRISTOL BS1 2NT 

(inducting the postcode) 



Patents ADP number (if you know it) 



05996483001 



If you are declaring priority from one or more 
earlier patent applications, give the country 
and the date of filing of the or of each of these 
earlier applications and (if you know it) the or 
each application number 



Country Priority application number 

(if you know it) 



Date of filing 
(day /month /year) 



If this application is divided or otherwise 
derived from an earlier UK application, 
give the number and the filing date of 
the earlier application 



Number of earlier application 



Date of filing 
(day /month /year) 



Is a statement of inventorship and of right 
to grant of a patent required in support of 
this request? (Answer Yes' if 

a) any applicant named in part 3 is not an inventor, or 

b) there is an inventor who is not named as an 
applicant, or 

c) any named applicant is a corporate body 
See note (d)) 



Yes 



Patents Form 1/77 



9. Enter the number of sheets for an^^Bb< 



following items you are filing with^lSPtbrm. 

Do not count copies of the same document 




l/Onunuauon sneeis ui uuo iuhu 
Description 


12 


Claims 
Abstract 


> 


Drawing^ 


- 


10. If you are also filing any of the following, 
state how many against each item. 

Priori tv documents 




Translations of priority documents 


- 


Statement of inventorship and right 
to grant of a patent (Patents Form 7/77) 


Two 


Request for preliminary examination 
and search (Patents Form 9/77) 


One 


Request for substantive examination 

(Patents Form 10/77) 




Any other documents 

(please specify) 


- 


11. 


I/We request the grant of a patent on the basis of this application. 

Signature ? <&J()^^ 07/02/03 
PAGE HARGRAVE 


12. Name and daytime telephone number of 
person to contact in the United Kingdom 


Mr WIJ Newstead (0117) 927 6634 



AlteSl application for a patent has been Sled, the Comptroller of the Patent OfBce will consider whether publication 
or communication of the Invention should be prohibited or restricted under Section 22 of the Patents Act 1977 You 
wOl be informed if it is necessary to prohibit or restrict your invention in this way. Furthermore if you hvein the 
United Kingdom, Section 23 of the Patents Act 1977 stops you hrom applying for a patent abroad without Brst getting 
written permission from the Patent OfBce unless an application has been Bled at least 6 weeks beforehand m the 
United Kingdom for a patent for the same invention and either no direction prohibiting publication or 
communication has been given, or any such direction has been revoked. 

a° ^Ifyou need help to Sll in this form or you have any questions, please contact the Patent OfBce on 0645 500505. 

b) Write your answers in capital letters using black ink or you may type them 

c) If there is not enough space for all the relevant details on any part of this form, please continue on a separate 
sheet of paper and write "see continuation sheet" in the relevant part(s). Any continuation sheet should be 
attached to this form. 

d) If you have answered "Yes' Patents Form 7/77 will need to be Bled. 

e) Once you have Blled in the form you must remember to sign and date it 

f) For details of the fee and ways to pay please contact the Patent OfBce. 



Patents Form 1/77 



DUPLICATE 



ESTIMATING THE TIME OF ARRIVAL OF A SEISMIC WAVE 

The present invention relates to a method of calculating an estimated tim e of arrival of a 
seismic or microseismicP or S wave at a sensor station, data processing means for carrying 
out such a method, and a data carrier comprising computer readable program means for 
carrying out such a method when run on a computer. 

A typical seismic or microseismic event produces both P-wave and S-wave signals. These 
distinct signals travel at different velocities through the earth. For a wide range of crustal 
rock types the ratio between the P-wave and S-wave velocity is well known and does not 
vary to any great extent. Indeed in many instances it is not possible a priori to obtain an 
accurate S-wave velocity and so the S-wave velocity structure for a region of interest is 
often derived by assuming axonstant ratio with the P-wave velocity structure. Thus we 
find that, in practice, the two seismic velocities structures within a volume of interest will 
show an approximately constant ratio to each other, even where the value of both types of 
velocity changes. 

The normal method of estimating the location and time (in terms of the coordinates x, y, 
z and time zero (tO)) at which a seismic or microseismic event occurs involves positioning 
a number of sensor stations, comprising geophones or the like, at various locations in a 
catchment -area, and measuring the arrival times at each of said stations of the P and S 
waves generated. Once the time of arrival of each wave at each station has been picked 
(i.e. decided upon based on the seismic data recorded at the station concerned), a guess is 
made at the location and time of the seismic or microseismic event, and the known or 
assumed P and S wave velocities are used to calculate the times at which the P and S . 
waves would reach each station, assuming that this guess were correct These calculated 
arrival times, based on the guessed location and time of origin, are then compared with the 
picked arrival times (which, as noted above, are based on the seismic data actually 
recorded) and the location and time of origin adjusted, via a mathematical processes known 
as the least squares process, such that the final estimated location and time of origin of the 
event best fits, in a least squares sense, with the picked arrival times of the waves at the 



stations. 



While the above method is generally accepted to be the best method of event location 
currently available, the reliability and accuracy of the location of a given seismic or 
microseismic event, recorded on a given network, is heavily dependent upon the accuracy 
and reliability of the arrival time picking. In some cases, the quality of the seismic data 
may make it difficult to identify a distinct arrival time at all. If a picked arrival time is 
badly incorrect, or worse still a wave type wrongly identified (for example where aP-wave 
picked as an S-wave or a P to S conversion picked as an S-wave), the final estimated 
location and time of origin of the event is likely to differ significantly from the true 
location and time. Furthermore, the least squares location process means that the effect of 
a single bad timing is spread out and manifests itself in all the channel residuals and the 
overall misfit of the location, making the identification of incorrect pick less than 
straightforward. 

According to one aspect of the present invention, data processing means for calculating an 
estimated time of arrival of a seismic or microseismic P or S wave at a sensor station are 
provided, comprising a data processor adapted to: 

a) calculate an estimated time of origin for. the seismic or microseismic event 
generating the P and S waves, based on a P to S wave velocity ratio and picked 
arrival times of the P and S waves at a sensor station other than the one for which 
the estimated time of arrival of the P or S wave is to be calculated; and 

b) calculate the estimated time of arrival of the P or S wave, based on a P to S wave 
velocity ratio, the estimated time of origin of the event and, where the estimated 
arrival time of a P wave is to be calculated, a picked arrival time of the S wave at 
the sensor station for which the estimated arrival time of the P wave is being 

, calculated or, where the estimated arrival time of a S wave is to be calculated, a 
picked arrival time of the P wave at the sensor station for which the estimated 
arrival time of the S wave is to be calculated. 



The data processor may be adapted to calculate estimated arrival times for both the P and 
S waves at a sensor station. 

Preferably, the data processor is adapted to calculate a plurality of estimated times of 
arrival of the P and/or S wave at a sensor station, based on a plurality of estimated times 
of origin for the microseismic event calculated from the picked arrival times of the P and 
S waves at a plurality of sensor stations other than the one at which the estimated times of 
arrival are to be calculated. The data processor may be further adapted to display the 
picked arrival times and estimated arrival times in relation to each other such that the 
clustering pattern of the arrival times can be analysed. Where this is the case, the data 
processor may also be adapted to display information regarding the calculation of any 
particular estimated arrival time in response to the selection of said estimated arrival time 
by an user. 

Preferably, the data processor is adapted to calculate one or more estimated times of arrival 
for the P and/or S waves at each sensor station in a network of sensor stations, wherein for 
each sensor station the necessary estimated time or times of origin are calculated from the 
picked arrival times of the P and S waves at one or more of the other stations in said 
network. 

In some cases, the data processor may require the user to pick arrival times for the P and 
S waves at the various sensor stations by studying the seismic data recorded at each station, 
and to then provide the data processor with the picked arrival times. Alternatively, the data 
processor may be adapted to receive seismic data 'from the sensor stations and to pick 
arrival times for the P and S wave at each sensor station based on said seismic data. 

Where a number of possible arrival times for a wave at a sensor station could be picked 
(either by the user, or by the data processor itself), the data processor may adapted to 
compare said possible arrival times with any estimates calculated for the arrival time of 
said wave at said station in order to determine which of the possible picked arrival times " 
are more likely to correspond to the true arrival time of said wave at said sation. The data 



processor may be further adapted to select or modify one of said possible arrival times in 
order to arrive at a final picked arrival time that, based on the above determination, seems 
most like to correspond to the true arrival time of said wave. Alternatively, the data 
processor may be adapted to indicate which of the possible arrival times should be selected 
or modified in order to arrive at a final .picked arrival time that seems, based on the above 
determination, to be most likely to correspond with the true arrival time. 

According to another aspect of the present invention, a data carrier comprising computer 
readable program means for adapting a computer to function as the data processing means 
according to the above aspect of Hie invention is provided. 

According to a further aspect of the present invention, a method of calculating an estimated 
time of arrival of a seismic or microseismic P or S wave at a sensor station is provided, 
said method comprising the steps of: 

a) calculating an estimated time of origin for the seismic or microseismic event 
generating the P and S waves, based on a P to S wave velocity ratio and picked 
arrival times of the P and S waves at a sensor station other than the one for which 
the estimated time of arrival of the P or S wave is to be calculated; and 

b) calculating the estimated time of arrival of the P or S wave, based on a P to S wave 
velocity ratio, the estimated time of origin of the event and, where the estimated 
arrival time of a P wave is to be calculated, a picked arrival time of the S wave at 
the sensor station for which the estimated arrival time of the P wave is being 
calculated or, where the estimated arrival time of a S wave is to be calculated, a 
picked arrival time of the P wave at the sensor station for which the estimated 
arrival time of the S wave is to be calculated. 

Thus the present invention provides means for obtaining one or more estimates as to the 
arrival time of the P and/or S waves at one or more of the sensing stations, without the 
location of the seismic or microseismic event first having to be located. The estimated 
arrival times thus produced can be used in a number of ways. Where the seismic data 



recorded at a station cannot be interpreted with sufficient reliability in order to make a 
pick, an estimate of the arrival time can be used to replace the pick that is not being made. 
Where a number of sensing stations are present, a spread of estimated arrival times can be 
calculated, in order to provide the user with sufficient information to identify possible bad 
picks before an attempt at event localisation is made. Such' spreads of estimates can also 
be used for analysis of the P to S wave velocity structure and ratio, as part of an auto- 
location process and/or for estimating picking uncertainty. 

A specific embodiment of the invention will now b e described, solely by way of example. 

According to the present embodiment, the data processing means comprise a computer 
running software for calculating the necessary estimated times of origin and arrival (though 
this could equally be achieved via suitable electronic hardware) by applying certain 
algorithms, the derivation and application of which are explained below. 

If one assumes that the network of sensor stations (each of which could be a seismic 
detector as described in GB-A-2 275 337 for example) are deployed in an arbitrarily 
complex isotropic earth with a constant or near to constant P- to S-wave ratio, then the 
distance (D) of a microseismic event from a station can be expressed as: 

D = (Ts - To) Vs 

and/or 

D = (Tp-To)Vp 

where Vs is the S-wave velocity, Vp is the P-wave velocity, Tp is the arrival time of the 
P-wave at the station, Ts is the arrival time of the S-wave at the station, and To is the time 
of origin of the microseismic event generating said P- and S-waves. 

Multiplying through the two equations gives: 



TsVs-ToVs = D 



and/or 

5 • TpVp-ToVp = D 

Eliminating D, the source to the distance, we get one equation: 
TsVs - ToVs = TpVp - ToVp 

10 

' Collecting all terms including To, the origin time, gives: 
To Vp - ToVs = Tp Vp - TsVs 
15 ••• ToVp/Vs - To = TpVp/Vs - Ts 
To(Vp/Vs - 1) = TpVp/Vs - Ts 

Vp/Vs represents the P-wave to S-wave velocity ratio, which we will designate as R. 
20 Substituting R for Vp/Vs gives: 

To®- l) = TpR-Ts 

To = (TpR-Ts)/(R-l) 

25 

We now have an estimate of the origin time of the earthquake in terms of the arrival times 
(Tp and Ts) of the P- and S- waves at the station and the ratio R between the P-wave to S- 
wave velocity field ratio. 

30 This equation has R in it twice and in order to re-arrange it so that R only appears once we 
can add a +Tp and a -Tp term to the top line of the equation, which gives: 



Re-arranging we get: 



To - 



Tp(R-l) + Tp - Ts 
(R-l) 



. Separating tenns gives: 



To 



Tp(R-l) + Tp _ _Ts 



(R-l) (R-D (R-l) 



Now we can cancel the (R-l) term in the first part of the right-hand side of the equation to 
give: 



To finally give: 

To = Tp-(Ts-Tp)/(R-l) 

Thus, where the arrival times of the P- and S-waves at a station can be picked with a 
reasonable degree of certainty and given a constant velocity-ratio earth mo del we can, from 
that station, get an estimate of the origin time of the event based solely on the picked 
arrival times of the P- and S- waves and the assumed P- and S- wave velocity field ratio. 
Note that this equations does not require knowledge of the actual location of the 
microseismic event or the velocity fields for the P- and S-waves. 

Where the station for which the P- and S- wave arrivals have been picked is designated as 



To = Tp + 



Tp Ts 



(R-l) (R-l) 




8 

station 1, we can add a 1 subscript to show that these terms apply to station 1 to get: 

To,-Tp,-Crs,-Tpi)/CR-l) C 1 ) 
5 where To , is the origin time estimated by using timings from station 1 . 

Likewise, if the arrival times were picked at a different station, designated station 2, we can 
add a 2 subscript to get: 

10 To^Tfe-CT^-Tp^CR-l) (2) 

Note that, given the assumption of an isotropic earth, the velocity field ratio remains 
constant for both stations. 

1 5 However, the times of origin calculated from stations 1 and 2 should be identical, since the 
same event is generating the waves recorded at both stations. Thus, in circumstances 
where an S-wave pick at station 2 has not been possible, or the accuracy of the pick seems 
uncertain, such that an estimate for the S-wave arrival time at station 2 is desired, this can 
be obtained by substituting To 2 for To, and Ts^ for.Tsz in equation (2) to get: 



20 



25 



30 



To^Tpj-CTs^-Tp^/CR-l) 

To,(R-l) = Tp 2 (R-l) - (Ts a - T P2 ) 
= Tp 2 (R-l)-T Se2 +Tp 2 

T Se2 = Tp 2 (R-l) - To,(R-l) +Tp 2 

= Tp 2 R - Tp 2 - To,R - To, + Tp 2 
= Tp 2 R-To,R-To l 
= <Tp 2 -To,)R-To, 

Ts e2 = R(Tp 2 -To 1 )-To 1 

•Where Ts^ represents an estimate for the S-wave arrival time at station 2. Note that the 



calculation of the estimated time of arrival of the S-wave at station 2 requires only the P 
and S picks from station 1 (in order to calculate To^, the P-wave pick at sation 2 and the 
assumed constant R ratio. 

Likewise, where an estimate for the P-wave airival time at station 2 is desired, this can be 
obtained by substituting To 2 for To,'and Tp^ for Tp 2 in equation (2) to get: 

To^lTp^-CTsa-TpJ/CR-l) 

T 0l (R-l) = T Pc2 (R-l) - (Ts 2 - T Pe2 ) 
= Tp e2 (R«l)-Ts 2 + Tp^ 

T 0l (R-l) + TS2 = Tp e2 (R-l) + T Pc2 

^Tp^R-Tp^ + Tp., 
= Tp e2 R 

Tp^^CTo^-^+Ts^/R 
^(To^-Toj -KTs^/R 
= TOjR/R + (Ts 2 - TOj)/R 

' Tp e2 = To, + (Ts i -To 1 )/R 
where Tp^ is the estimated P-wave arrival time for station 2. 

It follows, that where P-wave and S-wave picks for station 2 have been made, the process 
described above can be reversed in order to obtain estimates for the picks, both P-wave and 
S-wave, on station 1. Indeed, using the above process, every station for which a P-wave 
and an S-wave pick has been made allows an estimated time of arrival to b'e calculated for 
the P- and/or S-waves at another station in the network (assuming an arrival time for the 
P- or S-wave not being estimated has been picked for said station). 

Thus, for a network of six stations, if all the P-wave and S-wave picks are made at each 
station then five estimated times of arrival for each wave at each station can also be 
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calculated. If apick cannot be made for the arrival time of one wave at one of the stations 
(due to the arrival time of the wave not being identifiable from the seismic data) then one 
of the five estimated arrival times for said wave at said station (or an average thereof) can 
be used for the purpose of subsequent localisation of the microseismic event. Likewise, 
if all the picks have been made, but one of the picks is significantly different from the true 
arrival time (again due to the difficult of identifying the true arrival time from the data 
recorded), then the estimates produced from Hie bad pick will differ significantly from the 
other estimates which do not include the bad pick, allowing the bad pick to be identified 
and corrected (or indeed replaced by an estimate) prior to the localisation process. This 
is in contrast to attempting to examine the channel residuals produced by the localisation 
process in order to identify a bad pick, since each channel residue is influenced by all the 
picks used in the location as well as the P to S velocity structure ratio and three timings. 

The estimated arrival times that have been generated can also be used to help improve the 
consistency of arrival time picking. If the assumption of a constant P to S velocity ratio 
is true then all the estimated arrival times for a wave should form a tight cluster around the 
picked arrival time for a wave. If such a clustering pattern is not observed, there are three 
basic reasons why this may be the case, which are as follows: 

1) the assumption of a constant p to s ratio is . does not hold true for the area 
concerned; 

2) the p to s ratio used is incorrect; or 

3) one or moire of the picked arrival times is incorrect to a significant degree. 

Reasons (1) and (2) are related to the velocity structure. If the reason is reason (1) then this 
can be easily remedied by changing the P or S wave velocity until good clustering is 
observed. If the reason is reason (2) then the P to S ratio can be adjusted until Hie optimum 
ratio is found that gives the best clustering around a group of many timings from many 
events. Thus the process of estimating arrival times described herein could also be applied 
in order to optimise the velocity structure. Typically it is the S-wave velocity that is most 
poorly resolved and so in most cases it will be this velocity that needs to be altered in order 
to obtain the optimum velocity ratio and structure. It is possible to distinguish between 
reasons (1) and (2) in that where reason (2) applies (ie. where the ratio is constant but 
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inappropriate) good clustering together of the estimated picks will occur, but not be around 
the picked arrival times. Where reason (1) applies the effect will be more subtle and the 
estimated picks will not cluster to the same extent. In practice this may be a more difficult 
situation to identify, but in most cases the variations in the P to S ratio within a particular 
5 region are likely to be small, such that this will at most be a second order effect. Once 
many events have been located, itwouldbepossibleto investigate the relationship between 
. location and P to S ratio based on the data obtained. 

In most cases, reason (3) will be the cause of the desired clustering pattern not being 
1 0 observed. In order to assist the user in identifying which stations picked arrival times are 
giving the bad estimate, the software may further allow for the user to select any of the 
displayed estimates using the computer mouse or curser keys, in response to which the 
station that gave rise to the estimate is highlighted, hi this way the program implement on 
the computer guides the user to the station with the "worst" picks. This process can 
1 5 furthermore be repeated until the clustering of the estimated and measured pities appears 
satisfactory to the user. 

The program tunning on the computer may also be designed to use the estimated picks as 

part of an auto-location process. As has been described above, a "good" pick can be 

20 identified by the fact that it has a cluster of estimated picks around it. During the process 
_ 

of auto-timing it is common for the computer to initially select a number of potential picks 
for the arrival time of a wave at a station. By applying the above process of generating 
estimates for each pick, and analysing the clustering of the estimates, around each 
prospective auto-picks, the program can be designed to eliminate those prospective auto- 
25 picks that are less likely to correspond to the true arrival time. The program can also use 
the degree of clustering to decide if an event has been picked well enough to warrant 
storing and use for a subsequent localisation process. 

While estimating picking uncertainty or accuracy canbe a difficult and subj ective process, 
30 it is also possible to use the estimated arrival times in order to provide an indication of 
picking uncertainty. Each estimated arrive time includes two other arrival times. For 
example, the S-wave arrive time on station n involves the following picks: 
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Ts, 


T Pn 


1P2 


Ts, 


T Pn 


Tp 3 


Ts 3 


Tp n 


Tp 4 


Ts 4 


Tp n 


Tp 5 


Ts 5 


Tp„ 



La Which' Tpj indicates the picked arrival time for the P-wave at station 1, Ts x indicates the 
picked arrival time for the S-wave at station 1, Tp n indicates the picked arrival time of the 
P-wave at station n, and so on. From the above it will be apparent that Tp n is a constant 
10 in all the estimates. Thus the scatter in the estimates of the Ts n will be a function of the 
picking uncertainty in the P- and S- wave picking for each channel. 



CLAIMS 



1 . Data processing means for calculating an estimated time of arrival of a seismic or 
microseismic P or S wave at a sensor station, comprising a data processor adapted 
to: 

a) calculate an estimated time of origin for the seismic or microseismic event 
generating the P and S waves, based on a P to S wave velocity ratio and 
picked arrival times of the P and S waves at a sensor station other than the 
one for which the estimated time of arrival of the P or S wave is to be 
calculated; and 

b) calculate -the estimated time of arrival of the P or S wave, based on a P to 
S wave velocity ratio, the estimated time of origin of the event and, where 
the estimated arrival time of a P wave is to be calculated, a picked arrival 
time of the S wave at the sensor station for which the estimated arrival time 
of the P wave is being calculated or, where the estimated arrival time of a 
S wave is to be calculated, a picked arrival time of the P wave at the sensor 
station for which the estimated arrival time of the S wave is to be 
calculated. 

2. Data processing means according to claim 1 , said data processor being adapted to 
calculate estimated arrival times for both the P and S waves at a sensor station. 

3. Data processing means according to claim 1 or2, said data processor being adapted 
to calculate a plurality of estimated times of arrival of the P and/br S wave at a 
sensor station, based on a plurality of estimated times of origin for the 
microseismic event calculated from the picked arrival times of the P and S waves 
at a plurality of sensor stations other than the one at which the estimated times of 
arrival are to be calculated. 

4. Data processing means according to claim 3 , wherein the data processor is further 
adapted to display the picked arrival times and estimated arrival times in relation 
to each other such that the clustering pattern of the arrival times can be analysed. " 
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Data processing means according to claim 4, wherein the data processor is adapted 
to display information regarding the calculation of any particular estimated arrival 
time in response to the selection of said estimated airival time by an user. 

Data processing means according to any preceding claim, said data processor being 
adapted to calculate one or more estimated times of arrival for the P and/or S waves 
at each sensor station in a network of sensor stations, wherein for each sensor 
station the necessary estimated time or times of origin are calculated from the 
picked arrival times of the P and S waves at one or more of the other stations in 
said network. 

Data processing means according to any preceding claim, wherein the data 
processor is adapted to receive seismic data from the sensor stations and to pick 
arrival times for the P and S wave at each sensor station based on said seismic data. 

Data processing means according to any preceding claim, wherein, where anumber 
of possible arrival times for a wave at a sensor station could be picked, the data 
processor is adapted to compare said possible arrival times with any estimates 
calculated for the arrival time of said wave at said station in order to determine 
which of the possible picked arrival times are more likely to correspond to the true 
arrival time of said wave at said sation. 

Data processing means according to claim 8, wherein the data processor is further 
adapted to select or modify one of said possible arrival times in order to arrive at 
a final picked airival time that, based on said determination, seems most like to 
correspond to the true arrival time of said wave. 

Data processing means according to claim 8, wherein the data processor is adapted 
to indicate which of the possible arrival times should be selected or modified in 
order to arrive at a final picked airival time that seems, based on said 
determination, to be most likely to correspond with the true arrival time. 
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A data carrier comprising computer readable program means for adapting a 
computer to function as the data processing means of any one of claims 1 to 10. 



A method of calculating an estimated time of arrival of a seismic or microseismic 
P or S wave at a sensor station, said method comprising the steps of: 

a) calculating an estimated time of origin for the seismic or microseismic 
event generating the P and S waves, based on a P to S wave velocity ratio 
and picked arrival times of the P and S waves at a sensor station other than 
the one for which the estimated time of arrival of the P or S wave is to be 
calculated; and 

b) calculating the estimated time of arrival of the P or S wave, based on a P 
to S wave velocity ratio, the estimated time of origin of the event and, 
where the estimated arrival time of a P wave is to be calculated, a picked 
arrival time of the S wave at the sensor station for which the estimated 
arrival time of the P wave is being calculated or, where the estimated 
arrival time of a S wave is to be calculated, a picked arrival time of the P 
wave at the sensor station for which the estimated arrival time of the S 
wave is to be calculated. 
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ABSTRACT 



Data processing means are described for calculating an estimated time of arrival of a 
seismic or microseismic P or S wave at a sensor station, based on a P to S wave velocity 
ratio, a calculated estimated time of origin of the event and, where the estimated arrival 
time of a P wave is to be calculated, a picked arrival time of the S wave at the sensor 
station or, where the estimated arrival time of a S wave is to be calculated, a picked arrival 
time of the P wave at the sensor station. 



This Page is Inserted by IFW Indexing and Scanning 
Operations and is not part of the Official Record 



Defective images within this document are accurate representations of the original 
documents submitted by the applicant. 

Defects in the images include but are not limited to the items checked: 

□ BLACK BORDERS 

□ IMAGE CUT OFF AT TOP, BOTTOM OR SIDES 

□ FADED TEXT OR DRAWING 



□ SKEWED/SLANTED IMAGES 

□ COLOR OR BLACK AND WHITE PHOTOGRAPHS 

□ GRAY SCALE DOCUMENTS 

□ LINES OR MARKS ON ORIGINAL DOCUMENT 

□ REFERENCE(S) OR EXHIBIT(S) SUBMITTED ARE POOR QUALITY 

□ OTHER: 

IMAGES ARE BEST AVAILABLE COPY. 
As rescanning these documents will not correct the image 
problems checked, please do not report these problems to 
the IFW Image Problem Mailbox. 



BEST AVAILABLE IMAGES 




BLURRED OR ILLEGIBLE TEXT OR DRAWING 



